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ABSTRACT 

Many extended extragalactic radio sources require a local in situ acceleration mech- 
anism for electrons, in part because the synchrotron lifetimes are shorter than the bulk 
travel time across the emitting regions. If the magnetic field in these sources is localized 
in flux tubes, reconnection may occur between regions of plasma /3 (ratio of particle to 
magnetic pressure) << 1, even though jj averaged over the plasma volume may be > 1. 
Reconnection in low (3 regions is most favorable to acceleration from reconnection shocks. 
The reconnection X-point regions may provide the injection electrons for their subsequent 
non-thermal shock acceleration to distributions reasonably consistent with observed spec- 
tra. Flux tube reconnection might therefore be able to provide in situ acceleration required 
by large scale jets and lobes. 



Subject Headings : Magnetic Fields: MHD; Acceleration of Particles; Galaxies: Active, 
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I. Introduction 

The non-thermal emission from extended extragalactic radio sources is hkely syn- 
chrotron radiation by relativistic electrons (Begleman et al., 1984). Observationally, the 
flux satisfies F{iy) oc where a is the spectral index, and u is the frequency. This 

corresponds to a power law electron distribution with a number density N('ye) oc 7g"*, 
where 7e is the electron Lorentz factor and s = 2a + 1. For over 90% of the sources, 
0.5 < a < 1 (Begelman et al., 1984), corresponding to 2 < s < 3. Jets and hot spots show 
flatter spectra than diffuse radio lobes. 

Many large scale synchrotron sources have sizes larger than the distance electrons 
can be convected in a synchrotron loss time, exhibiting the need for an in situ accelera- 
tion process (Achterberg, 1987). When the magnetic field is sufficiently tangled, electron 
transport is due to the bulk flow. The transport distance in a synchrotron loss time is 

Dsyn = V.Tsyn = V^^-i^m^C^ / P^yn = (8 X 10«)H/7eS^ (1) 

where Psyn is the synchrotron power, rUe is the electron mass, B is the magnetic fleld 
magnitude, and Vj, is the bulk flow speed. An electron of energy m^.c^je radiates according 
to (Rybicki & Lightman, 1979) 

u^l.Gx 10-2(S/10-^G)(mec27e/GeV)2 GHz. (2) 

Using (1) and (2), an electron in a Mpc scale lobe radiating at 1.5 GHz in a 5 /xG fleld 
has 7e ~ 7.7 x 10"^, and will lose its energy in > 10 Mpc if Vf, = c. However, if a 
thermal plasma is present in the lobes, plasma instabilities can reduce the bulk veloc- 
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ity to V5 < c (Achterberg, 1987). For (3 = P/Pb < 1, where P is the particle pres- 
sure and Pb is the magnetic field pressure, the maximum becomes the Alfven speed 
Va ~ 5 X 10^(S/5 X 10~^G)(n/10~^cm~^)~-^/^ cm s~-^, where n is the number density. 
Thus, Dgyn < 20 kpc. Even more suggestive are observations of optical synchrotron jet 
emission, for example in 3C33 (Meisenheimer and Roser, 1986), which imply ~ 100 GeV 
electrons (7 ~ 10^) radiating over several kpc in an equipartition (EQP) field of 10~^G. 
The synchrotron lifetime is then ~ 200yr, much shorter than the jet length travel time, 
so that in situ acceleration is required for the jet emission to last. HST observations sup- 
port this (Macchetto, 1992). In addition, many jets do not show the spectral line index 
gradients expected from expansion or synchrotron losses (Achterberg, 1987). 

The energy for particle acceleration and turbulence is converted from the bulk fiow. 
Because the turbulent motions are responsible for stretching the magnetic field, the average 
magnetic field energy will be at most equal to the average turbulent energy. These points 
are summarized for primarily non-relativistic bulk fiows by 

(l/2)p«„eH' > {imPaveVl > {B^)/87T and {l/2)paveVb^ > P ~ Pec\le), (3) 

where pave ^tnd pe ^tre the average and electron densities, Vl is the outer eddy scale velocity 
and the brackets indicate the average. The presence of a large scale mean fiow and mean 
magnetic field does not preclude the presence of a tangled small scale field. In the Galaxy 
for example, the rotation provides an underlying ordered fiow, while the observed large 
and small scale fields have the same magnitude (Heiles, 1995). 

Previous studies have appealed to stochastic and fast shock acceleration (Achterberg, 
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1987; Eilek & Hughes, 1993). The first results directly from turbulent motions in magne- 
tized plasmas or from low frequency plasma and MHD waves. The second occurs near fast 
shocks that may be present in the flow. It is possible for both mechanisms to accelerate 
particles to power law spectra. However, both are "re-acceleration" mechanisms; electrons 
need to be pre-accelerated before these processes can maintain high energy electron tails. 

In this paper, a kind of hybrid approach is suggested which appeals to slow shocks that 
are naturally produced from magnetic reconnection events (for an alternative approach to 
reconnection in jets, see Romanova & Lovelace, 1992). The shocks would occur throughout 
the plasma, wherever a reconnection site occurs. Reconnection shocks also have the unique 
feature of an adjacent X-point region where the field annihilates. This region may be able 
to directly inject electrons to their required pre-accelerated energies. Section II addresses 
how reconnection shocks can accelerate electrons. The most efficient shock acceleration, 
with power law indices appropriate for large scale jets, would occur when the reconnection 
takes place in regions of low plasma (3 (Blackman & Field, 1994). In section III, I point 
out that this may occur in an astrophysical plasma even if the average ratio of the particle 
to magnetic pressure Pave ^ Ij when the reconnection occurs between reduced density fliix 
tubes. In Section IV, I address the application to jets and lobes more specifically, and 
summarize in section V. 

II. Reconnection Shock Acceleration 

Magnetic reconnection occurs as regions of skewed magnetic polarity intersect. The 
intersection produces a thin dissipation region where flux freezing is violated and magnetic 
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field is annihilated (e.g. Biskamp, 1994). The annihilation produces a topology change 
with an "X-point" at the interface. Much of the work in reconnection theory has focused 
on determining how fast oppositely magnetized plasma flows can merge across an X-point 
region in the steady-state. Although the dependence of this rate on Rm has been debated, 
many simulations show the presence of slow mode shocks (Biskamp, 1994), highlighting 
the abrupt change in direction between the inflow to the X-point and outflow from the 
X-point. Such slow shocks are distinguished from fast mode shocks in part by the fact 
that the magnetic field downstream of a slow shock has a lower magnitude than its value 
upstream — the opposite to the fast shock. 

Blackman & Field (1994), show that low (3 reconnection slow shocks are the strongest 
slow shocks, giving compression ratios q of 2.5 < q = pd/pu < 4, where Pd{u) is the 
downstream (upstream) density (see also Kantrowitz & Petschek, 1966). The upper limit 
corresponds to a quasi-parallel shock (upstream magnetic field nearly parallel to the shock 
normal), and the lower limit corresponds to a quasi-perpendicular shock (upstream mag- 
netic field nearly perpendicular to the shock normal). The importance of q can be seen 
from the analytical approach to shock acceleration. Assuming that gradients in the nor- 
mal direction >> those along the shock, the diffusion-convection equation across a shock 
is given by (Jones & Ellison, 1991) 

dn[VnN - kAN] - {l/3){dMdp^\ppN] = 0, (4) 

where Vn is the normal flow velocity across the shock, Kn is the normal diffusion coefficient, 
and Pp is the particle momentum. Fermi acceleration operates as the particles diffuse 
between scattering centers (presumably Alfven turbulence) on each side of the shock. 
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Particles always see the centers converging, as the normal velocity is larger upstream. The 
solution of (4) across the shock with thickness << mean free path (Jones & Ellison, 1991) 
indeed shows that the outflow energy spectrum for a steeper inflow spectrum takes the 
power law form N oc 7"*, with energy index s = (g + 2)/(g — 1). For 2.5 < q < 4 we 
then have 2 < s < 3 for these shocks, which is consistent with the observed range for 
extragalactic radio sources given in section I. 

Though an analytical treatment was employed above, shock acceleration is actually 
a very non-linear process; the Fermi acceleration engine across shocks is very efficient, 
transferring > 1/5 of the inflow energy to particles (Jones & Ellison, 1991). These particles 
tend to smooth out the shock by diffusion, produce turbulence, and can even increase 
the compression ratio above the jump condition value, as their escape and acceleration 
change the downstream equation of state. Non-thermal acceleration is enhanced in the 
non-linear regime. For an ion-electron plasma, the signature of effective non-linear shock- 
Fermi acceleration in simulations is the appearance of electromagnetic beam instabilities 
(Jones & Ellison, 1991), brought on by the interaction of back-scattered, energized, ions 
with the inflowing plasma. The instabilities indicate that non-thermal ions are produced 
from an initially thermal input, initiating the Fermi acceleration process. Such instabilities 
have been seen in fast as well as in slow shock simulations (Omidi & Winske, 1994). 

Some observational support for reconnection shock acceleration is present in the geo- 
magnetic tail where turbulence, required for Fermi acceleration, is seen on both sides of the 
shock fronts (Coroniti, et al., 1994). In addition, although non-thermal ions are observed. 
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non-thermal tails in the electron spectra are also seen (Feldman et al., 1990). These have 
not been modeled by hybrid simulations which assume a fluid electron population. 

It is indeed the acceleration of electrons, not ions, that is of interest for radio sources. 
Simulations show electron acceleration, but only for electrons that are injected above some 
critical energy (Ellison, 1992). Slow shocks would operate similarly with respect to Fermi 
acceleration. (For relativistic flows, Fermi acceleration would only operate for parallel 
shocks, while for perpendicular shocks, non-thermal electron acceleration requires a small 
proton fraction even for a primarily electron-positron jet [Hoshino et al., 1992]). Injection 
is needed because Fermi acceleration demands that the downstream particles be able to 
diffuse upstream. Particles must have enough energy to resonantly interact with the plasma 
waves that provide pitch-angle scattering. For Alfven turbulence (Eilek & Hughes, 1991), 
the electron lower bound is a factor ~ mp/me times that for protons and is given by 

7e > 1 + (mp/me) iVAdlcf (5) 

for < {me/mp)'^c'^ while for > {me/mp)^c^ , 

7e > {mp/me){VAd/c), (6) 

where rup is the proton mass, and VAd is the downstream Alfven speed. 

For reconnection shocks, the associated highly dissipative energy conversion near the 
X-point may provide the injection electrons. To see this, note that upon absorbing anni- 
hilated held energy at the X-point, the average 7e there ~ 1 -|- (y|^/c^)(mp/me), where 
Vau is the upstream Alfven speed. (Note that when B^^/Airp > c, then Vau is replaced by 
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IAvYau where ^au is the associated Lorentz factor.) Since ^avVau > lAdVAd for ^ slow 
shock (Blackman & Field, 1994), we see that an X-point region can therefore in principle 
always inject for both cases (5) and (6). 

To solve the in situ acceleration problem for radio sources, shock acceleration must 
dominate synchrotron loss. Thus, Tayn must exceed the shock acceleration time Tgh'- 

Tsyn = lefrieC^/Psyn = 67rmeC/(7eS^cr) > Tsh ^ ^n/V^, (7) 

where Vy, is the upstream velocity, a is the Thomson cross section, and is the diffusion 
coeflBcient normal to the slow shock. For electrons moving at c, k„ ~ c(7emeC^/eSd) 
(Achterberg, 1987). From (7), the condition for Tgyn > Tgh (also justifying the absence of 
a synchrotron loss term in (4)) is then 

7e < 3.8 X 10^(K/10^ cm/s)(S„/10-^Gauss)-^/2Cos6' (8) 

where ~ BuCos9 (Blackman & Field, 1994), has been used, and 9 is the angle of the 
upstream field with respect to the shock normal. 

III. Role of Flux Tubes for /3ave > 1 

Shock acceleration as described above, may occur between regions of low /? plasma 
even when Pave ^ 1; if the field is confined to fiiix tubes. Vishniac (1995) considers 
the steady-state dynamics of flux tubes in a turbulent plasma with (1/2) pV^ ~ (S^)/87r 
and Pave >> 1- In this case, the flux tubes flU a negligible fraction of the total volume. 
However, if the relations in (3) are near equalities, then (3ave ~ 1- Here the flux tubes 
would flU 1/2 the total volume. As this regime is outside of the scope of Vishniac (1995), 
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it will be considered in more detail below. It must be noted that in neither Vishniac 1995, 
nor the present paper, is the formation of flux tubes solved as an initial value problem. 
The presence of a steady-state is assumed and the consequences are explored. 

For jets, Kelvin-Helmholtz instabilities between the jet flow and the ambient medium 
(eg. Birkinshaw, 1991) are a source of turbulent eddies. Each energy containing (outer) 
scale eddy of wavelength L would stretch a tube to length L and radius r^, and the steady- 
state plasma would then contain a mesh of reconnecting fliix tubes. The thickness of each 
tube, r^, can be estimated by balancing the magnetic and turbulent eddy drag forces (to 
be justifled later) (Landau & Lifshitz, 1987; Vishniac, 1995)^. This gives 

(S2/47rre)(7rr2) ^ CdPextVl2ru (9) 

where Bt is the magnitude of the held in the fliix tube, pext is the density outside the fliix 
tube, and is the coefficient of turbulent drag. Since L is a wavelength, the radius of 
curvature Tc can be estimated by L/4 when the tube maximally responds to the turbulence. 
In equilibrium, i?|/87r ~ Pexti where P^xt is the external pressure, so that (9) gives 

n = LCdT Ml/ 477 ~ LCd/27r ~ 0.4L/27r ~ L/16, (10) 

where M| = Vl / {VP^^t/ Pext) ~ M| ~ 2/r, in EQP and V is the adiabatic index. For the 
last similarity in (10), Cd was estimated from the "drag" crisis (Landau & Lifshitz, 1987) 
which reduces < 1 at large Rl. Assuming > 1000, acceptable for jets and lobes 
(Begelman et al., 1984), Cd ~ 0.4. 

The steady-state MHD solution for the interior of a fliix tube with large aspect ratio 
(Vishniac, 1995) shows that the density in a tube falls off rapidly from the outside within 
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a thin shell of thickness ~ [i^MTedY^'^ ~ ^I^m ■> where Ted = L/Vl-, i^m is the magnetic 
diffusivity, and Rm is the magnetic Reynolds number. Inside the skin layer, the Alfven 
speed increases rapidly toward the center, due to the reduction of the tube central density. 
The tube reconnection time scale. Tree-, is the time to reconnect the skin layer plus the 
time to reconnect half the internal section of a tube. Since «rt we have 

Tree = F{RM){rs/VAs + n/VAi) ~ F {Rm)[L / {R^j!j^Vas) + L / {IWai)], (11) 

where Vai = Vau is Alfven speed inside the tube, Vas ~ Vl is the Alfven speed within r^, 
and F{Rm) is either for Sweet-Parker reconnection, or Logi?^ for Petschek recon- 
nection (Parker, 1979). We need to estimate Vai, and thus the mass fraction in the tubes, 

to obtain Tree and Kec = K = n/Trec- 

To estimate the tube reconnection velocity, and show that (3t « 1, justifying the use 
of the strong compression ratios as given in section II, it is first necessary to distinguish 
between the Pave » 1 case of Vishniac (1995) with the Pave ~ 1 case likely relevant 
for jets. In the former, tubes are nearly completely evacuated, and far from their nearest 
neighbor while for the latter, each flux tube would be in contact with ~ L/2rt ~ 8 other 
tubes because approximately 1/2 the volume is filled with fiux tubes. The magnetic field 
structure would therefore show a close mesh of reconnection sites. Each reconnection 
site occurs where a tube flattens up against a partner. The dissipation region at each 
site has a volume of Wrr ~ (2rs)(2rt)^. Because magnetic field is annihilated in each 
dissipation region, some material from the outside can leak in and slip onto tube field lines 
there; every field line in the tubes must eventually pass through X-point regions. A full 
reconnection event requires a height of 2rt to reconnect, so the amount of matter that 
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loads into a given tube during each Tree is ~ NrrPextWrri'^rt/2rs), where Nrr ~ 8 is the 
number of reconnection sites per tube. Since reconnection of closed loops removes 1/2 of 
the matter in any tube during Tree, only 1/2 of the loaded matter remains. The steady- 
state mass fraction of matter in the flux tubes is then {l/2)NrrPextWrr{'2rt/'2'fs)/[pextL^] ~ 
32r'l/L^ ~ 4x 10""^, using (10). Since the tubes comprise 1/2 the total volume, the density 
in the tubes is then ~ 4 x 10~^ pexti and though the tubes are not completely evacuated, 
Pt ^ 4: X iO~^Pave << 1) where Pt is P inside the tubes. Given the density reduction in 
the tubes, Vai = Vau ~ (4 x 1[)-^)-^/'^Vas ~ 16Vl, so from (11) 

Vrec = K ~ n/Trec ~ n/(TedF(i?M)/100) ~ (12) 

where the last equality follows for Petschek reconnection with F{Rm) ~ Log[i?M] ~ 25 
and using (10). The tube reconnection, which occurs between regions with Pt « 1, can 
in principle proceed with Vree ~ even for Pave ~ 1- That V^ec ~ supports the use 
of (11) even if tubes fill a large volume fraction. 

Note that flux tube mobility requires that the turbulent damping viscosity be neg- 
ligible (Vishniac 1995). This condition is vt < VLrt/n^ where ut is the viscosity. For 
Pave ~ 1: this requires Rl = VlL/ut > n^{L/rt) ~ 500. Rl ~ 1000 was used above. 

IV. Application to Extragalactic Sources 

Solving the in situ acceleration problem by the flux tube scheme requires that 1) (8) 
is satisfied and that 2) electrons encounter a reconnection acceleration site within their 
synchrotron loss times. Consider requirement 1: Since ~ Vree ^ Vl ^ Vb, using (12) 
and EQP, and B u ~ Bf, we can use (8) and (2) applied to flux tubes to obtain the new 
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condition 



1.6 X l{)'^{v/GHzY''^ < 10^(V6/10S cm/s)Cos^, 



(13) 



which is satisfied for reasonable choices of jet and lobe parameters. 

Requirement 2 means that the distance between sites must be less than Dsyn, or more 
conservatively, that every line of sight must pass through a reconnection region. Since 1/2 
of any tube of length L is incurring reconnection for Pave ~ 1 as described above, we have 
Dj > 2L ~ 8rt where Dj is the jet diameter. This provides an upper limit on the size of 
the turbulent outer scale given by L < Dj/2. If Dj ~ O.Skpc for example, this implies 
L < 250pc. From (6), this in turn implies that rt ~ 15.6pc. The size scale of tubes in the 
jets is a lower limit to the size scale of tubes in the lobes. 

Because L is much smaller than jet lengths of interest, a jet field would have many 
reversals along its length. This is consistent with observations in jets if EQP is assumed 
(Begelman et al., 1984). Field lines parallel to a jet can carry a magnetic fiiix < lO^^G 
cm^ near a central black hole of IO^Mq and then < lO^^G cm^ at Ipc. However, the 
observed EQP field strengths and cross sectional areas at > 10 kpc imply more like 10^°G 
cm^. Thus the field likely reverses many times, consistent with the fiux tube scheme. 

The total luminosity produced by the fiux tube model is given by 



where Lrr is the luminosity per reconnecting region. Near EQP, Lrr ~ 4x10 ^{B^/87r)Wrr{Vrec/rt)M 



Ltot — NrrL^ 



(14) 



12 



Then using (10) and (12), (14) becomes 

Ltot ~ 4 X 10-3(s2/87r)W^tot(10/Ted), (15) 

where Wtot is the total volume of the region of interest. For a lobe of diameter IMpc, Bt ~ 
5^G, and ~ lO^cm/s, Ted L/Vl-^ 250pc/H ~ 7.5 x lO^^sec so Ltot ~ 5 x lO^^erg/sec. 
For kpc scale jets with length 10 kpc, width 0.5 kpc, lOOGeV electrons, Bt ~ 10~^G, and 
Vi) ^ Vl ^ 0.1c, Ltot ~ lO'^'^erg/sec. These are within the range of observed luminosities 
for powerful radio sources (Muxlow & Garrington, 1991). 

V. Conclusions 

That slow shocks from reconnecting flux tubes might be able to provide the in situ 
acceleration of electrons required in large scale jets and radio lobes is suggested. If a steady- 
state flux tube structure ensues, reconnection regions will occur naturally throughout the 
turbulent plasmas, and therefore, so will the associated slow shocks. Even if (3ave ^ Ij 
the flux tubes are low /? regions. The slow shocks may then Fermi accelerate electrons 
to power law spectra reasonably consistent with the observed ranges for radio sources. 
The standard electron injection problem for shocks may be overcome by direct injection 
from the X-point regions. Future work should include a more dynamical approach, and a 
consideration of relativistic bulk flows. 

1) Eqn. (9) differs from Vishniac (1995) by a factor of 2. 
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